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ABSTRACT 

While  entanglement  is  believed  to  underlie  the  power  of  quantum  computation  and  communication,  it  is  not  generally 
well  understood  for  multipartite  systems.  Recently,  it  has  been  appreciated  that  there  exists  proper  no-signaling 
probability  distributions  derivable  from  operators  that  do  not  represent  valid  quantum  states.  Such  systems  exhibit 
supra-correlations  that  are  stronger  than  allowed  by  quantum  mechanics,  but  less  than  the  algebraically  allowed 
maximum  in  Bell-inequalities  (in  the  bipartite  case).  Some  of  these  probability  distributions  are  derivable  from  an 
entanglement  witness  W,  which  is  a  non-positive  Hermitian  operator  constructed  such  that  its  expectation  value  with  a 
separable  quantum  state  (positive  density  matrix)  psep  is  non-negative  (so  that  Tr[  W  p]<  0  indicates  entanglement  in 
quantum  state  p).  In  the  bipartite  case,  it  is  known  that  by  a  modification  of  the  local  no-signaling  measurements  by 
spacelike  separated  parties  A  and  B,  the  supra-correlations  exhibited  by  any  W  can  be  modeled  as  derivable  from  a 
physically  realizable  quantum  state  p.  However,  this  result  does  not  generalize  to  the  n-partite  case  for  ri> 2.  Supra- 
correlations  can  also  be  exhibited  in  2-  and  3 -qubit  systems  by  explicitly  constructing  “states”  O  (not  necessarily  positive 
quantum  states)  that  exhibit  PR  correlations  for  a  fixed,  but  arbitrary  number,  of  measurements  available  to  each  party. 
In  this  paper  we  examine  the  structure  of  “states”  that  exhibit  supra-correlations.  In  addition,  we  examine  the  affect  upon 
the  distribution  of  the  correlations  amongst  the  parties  involved  when  constraints  of  positivity  and  purity  are  imposed. 
We  investigate  circumstances  in  which  such  “states”  do  and  do  not  represent  valid  quantum  states. 

Keywords:  quantum  entanglement,  quantum  non-locality,  non-signaling  theories,  Popescu-Rohrlich  boxes,  EPR 


1.  INTRODUCTION 

Physics  imposes  limits  on  the  correlations  that  can  be  observed  by  distant  (i.e.  spacelike  separated)  parties.  In  particular, 
special  relativity  (SR)  implies  the  principle  of  no-signaling  (NS),  that  is  correlations  cannot  lead  to  any  sort  of 
instantaneous  communication  between  spacelike  separated  observers.  Quantum  correlations  may  be  stronger  than 
classical,  and  their  violation  of  Bell  inequalities1  (BI)  suggest  that  quantum  mechanics  (QM)  cannot  be  regarded  as  a 
local  realism  theory.  Tsirelson2  showed  that  there  is  an  upper  bound  to  the  violation  of  BI,  which  implies  that  the  amount 
of  non-locality  allowed  by  QM  is  limited.  Popescu  and  Rohrlich  (PR)  showed3  that  there  exists  a  broad  class  of  no¬ 
signaling  theories  which  allow  stronger-than-quantum  or  supra-quantum  correlations.  PR  developed  a  valid  joint 
probability  distribution  whose  violation  of  the  BI  lie  above  those  of  physical  quantum  correlations  and  below  the 
allowed  algebraic  maximum  of  the  BI  (the  latter  are  called  PR-Boxes).  Thus,  the  principle  of  NS  imposed  by  SR  does 
not  single  out  QM  from  these  other  post-quantum  NS  theories4  (PQNS). 

These  PQNS  have  much  in  common  with  QM  such  as  no-cloning,  information-disturbance  tradeoffs,  security  for  key 
distribution,  and  others.  Recently,  van  Dam5  showed  that  PR-Boxes  make  communication  complexity  trivial,  which  is 
not  the  case  within  QM.  Other  researchers  have  shown  that  PQNS  theories  would  lead  to  implausible  simplification  of 
distributed  computational  tasks  (see  Pawlowski6  and  references  therein).  It  is  now  widely  believed  that  theories  in  which 
communication/computational  complexity  is  trivial  are  very  unlikely  to  exist.  It  is  therefore  important  to  understand  the 
structure  of  the  PQNS  and  ultimately  to  find  physical  and  informational  principles  that  rule  them  out.  In  this  paper  we 
take  steps  in  that  direction  by  investigating  the  structure  of  PR  correlations  by  forming  operators  which  reproduce  these 
PR  probability  distributions.  We  investigate  circumstances  in  which  they  do  and  do  not  represent  valid  quantum  states. 
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2.  BELL  INEQUALITIES,  PR  BOXES  AND  SUPRA-QUANTUM  CORRELATIONS 

2.1  Bell  Inequalities  (BI) 

Nonlocality  is  expressed  by  means  of  violations  Bell  inequalities1  (BI)  which  set  upper  bounds  for  classical  correlations 
arising  from  local-realistic  theories.  For  bipartite  systems,  the  most  well  know  BI  is  the  Clauser-Horne-Shimony-Holt 
(CHSH)  inequality7  defined  as  follows.  Consider  a  bipartite  system  AA  ®2,  Alice  and  Bob,  each  possessing 
measurement  directions  A,B  =JA  and  C,D  =2  taking  measurement  values  a,b,c,d  ={±  1}.  We  define  the  correlation 
E(AC)  between  A=AA  and  C=2  as 


E(AC)  =  (AC)  =  £  acP(a,c\A,C) 

0,c={±  1}  (1) 

=  P(+,  +  I  A,C)  +  P(-  -  [  A,C)~  P{+,  -  |  A,  C)  -  P(-  +  I  A,C) 

In  (1),  we  define  P(a,c\A,C)  as  the  joint  probability  that  given  the  (inputs)  measurement  directions  A  for  Alice  and  C  for 
Bob,  Alice  obtains  the  (outputs)  measurement  result  a  and  Bob  obtains  the  value  b ,  subject  to  the  normalization 
condition  ^  _i+[ tP(a,c  \  A,  C)  =  1,  \/A,C.  Finally,  we  define  the  following  CHSH  correlation  parameter  S  by 

S  =  E(AC)  +  E(BC)  +  E(BD)-E(AD).  (2) 

S  has  been  cleverly  constructed  as  the  expectation  value  of  the  quantity  Arg  =A(C-D)  +  B(C+D).  If  A,B,  C,D  are  classical 
random  variables  taking  values  ±  1  then  it  can  be  readily  seen  that  if  (i)  C=D ,  then  |Arg|  =  \B(2C)\  =  2  and  if  (ii)  C=-D , 

| Arg|  =  \A(2D)\  =  2.  Thus,  for  classical  correlation  we  have  the  CHSH  inequality 

CHSH  inequality:  |  S  =  E(AC)  +  E(BC)  +  E(BD)  -  E(AD)  \<  Sa  =  2  (3) 

(where  the  subscript  “Cl”  denotes  “classical”).  For  a  large  class  of  measurement  directions  (but  not  all),  quantum  states 
can  violate  the  CHSH  inequality  (i.e.  |S|>2)  up  to  a  maximum  value  shown  by  Tsirelson2  to  be  Sq  =  2v/2.  Here,  a 
quantum  state  is  defined  as  a  positive  (i.e.  non-negative  eigenvalues)  Hermitian  matrix  with  unit  trace  denoted  by  the 
symbol  p.  The  archetypical  example  is  the  singlet  (Bell)  state 

Aingle,  =(|tl)-|lt))/V2  =(|01)-|10»/V2  (4) 

with  measurement  directions  in  the  x-y  plane:  A  =  x,B  =  y,C  =  (x  +  y)/s[2  ,D  =  (x-y)/s/ 2  that  saturates  the  Tsirelson 
bound  with  S=  -Sq  =  -2v/2.  This  is  a  manifestation  of  the  stronger  than  classical  correlations  that  can  be  exhibited  by 
quantum  states.  (Note:  quantum  states  with  measurement  directions  such  that  the  CHSH  inequality  is  satisfied,  i.e.  S< 2, 
are  not  distinguishable  from  classical  states  by  the  correlation  parameter  S ). 

It  is  instructive  to  note  that  the  CHSH  inequality  in  (3)  can  be  derived8  as  a  statement  of  a  classical  quadrilateral 
inequality  for  the  correlation  metric  A  (AC)  - 1  -  E(AC)  =  P(+,  -  |  A,  C)  +  P(-,  +  \  A,  C)  >  0.  Substituting  this  expression 
into  (3)  yields  A  (AC)  +  A  (BC)  +  A  (BD)  >  A  (AD)  =>  S  =  E(AC)  +  E{BC)  +  E(BD)  -  E{AD)  <  +2  (see  Fig.  1).  Thus,  the 

For  a  bipartite  system  AA  s 'B  with 
measurement  directions  A,B  =  A  and  C,D  £'B 

Classically:  A(AC)  +  A(BC)  +  A(BD)  >  A(AD) 

Quantum  Mechanically:  this  can  be 
violated, e.g.  A(AC)  +  A(BC)  +  A(BD)  <  A(AD) 
with  singlet  state  (l^)— l-i-T))/-^ =(|oi)— |io))/V2 

Fig.  1  CHSH  inequality  derived  as  a  violation  of  the  classical 


S  J4 

c  A(BC)  B 


quadrilateral  inequality. 
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violation  of  the  CHSH  inequality  by  quantum  states  can  be  interpreted  as  a  violation  of  the  classical  quadrilateral 
inequality  which,  for  certian  measurement  directions,  yields  the  distance  A  (AD)  via  the  direct  path  A-D  to  be  smaller 
than  the  sum  of  the  distances  around  the  indirect  path  A-C-B-D. 

Returning  to  the  CHSH  inequality  (3),  one  notes  that  it  is  bounded  by  the  algebraic  maximum  |S|  <  Sam= 4.  This  follows 
from  the  fact  that  the  correlations  E  are  bounded  by  |E|<1.  This  latter  result  can  be  inferred  by  writing  E  =  P++  +  P__  - 
(P+_  +  P_+)  =  2 (P++  +  P_)- 1  =  1-2 (P+_  +  P.+),  where  =  P++  +  P..  +  P+.  +  P.+  =1  has  been  used.  Using  the  fact  that  0<  P++ 
+  P-  <1  and  0<  P+_  +  P.+  <1  in  the  previous  two  expressions  for  E,  yields  the  desired  bound  |£]<1.  Therefore,  if  the  first 
three  correlations  in  (3)  take  the  value  ±1  and  the  last  correlation  takes  the  value  +1 ,  we  obtain  S=±4.  The  implication  of 
this  observation  is  that  the  regime  2a/2  <  S  <  4  represents  supra-correlations  that  are  stronger  than  quantum,  yet  are 
unphysical  by  Tsirelson’s  bound,  i.e.  cannot  be  realized  by  any  physical  quantum  state.  The  salient  question  to  study  is 
what  ‘natural’  principles  determine  the  exclusion  of  such  supra-correlations.  As  a  first  hypothesis,  one  might  surmise 
that  the  principle  of  no-signaling  from  special  relativity  (i.e.  that  information  cannot  be  instantaneously  broadcast 
between  spacelike  separated  observers)  might  exclude  supra-correlations.  Surprisingly,  this  is  not  the  case.  In  1994, 
Popescu  and  Rohrlich3  (PR)  were  able  to  construct  a  valid  joint  probability  distribution  between  a  pair  of  spacelike 
separated  observers  that  (i)  satisfies  the  non-signaling  principle,  and  (ii)  yields  the  algebraic  maximum  correlations 
allowed  by  the  CHSH  inequality.  Here  the  adjective  ‘valid’  implies  that  the  joint  probability  distribution,  and  all  its 
derived  marginal  probability  distributions  obtain  values  between  0  and  1 ,  and  satisfy  the  appropriated  normalization 
requirements  (i.e.  the  joint  and  all  marginal  probability  distributions  summed  over  all  outcomes  for  any  measurement 
settings  yields  unity).  These  correlations  are  now  called  PR  correlations,  which  we  describe  in  the  next  section. 

2.2  No  Signaling  (NS)  Theories  and  PR  Correlations 

We  wish  to  consider  correlations  between  n  spacelike  separated  parties  (observers)  Ah...^4„,  who  can  perform  m 
possible  measurements  xh...,x„  (x;={0,l,...,»!-l),  with  r  possible  outcomes  au...,an  (ai={0,l,...,r-l)  The  observed 
correlations  will  be  described  by  the  joint  probability  distribution  P(al,a2,..,,an  \x[,...,xn)  giving  the  probability  that 
the  parties  obtain  the  measurement  values  (outputs)  al,...,a„  when  their  local  measurement  apparatuses  (inputs)  are  set  to 
x\,...,x„.  The  joint  probability  distribution  is  constrained  only  by  the  conditions  0  <  P(al,a2,...,an  \  xl,...,xn)  <  1  and  the 

normalization  condition  V  P(a,,a,,...,a  \  x,,...,x)  =  1  for  all  measurement  settings Xi,...pc„. 

*-^a\ 

Imposing  the  no-signaling  (NS)  constraint,  i.e.  adherence  to  the  requirement  from  special  relativity  that  spacelike 
separated  measurements  should  not  influence  each  other  due  to  the  finite  speed  of  light  (communication),  requires  that 
the  marginal  probability  distributions  satisfy  the  additional  condition 

No  Signaling:  P(al,a2,...,ak  \  x1,...,xn)  =  ^  P{al,...,an  \  xx,...,xn)  =  P(ai,a2,...,ak  \  xx,...,xk).  (5) 

“m . “.ef0-1} 


Here,  the  first  equality  in  (5)  formally  defines  the  marginal  probability  distribution  describing  the  measurement 
outcomes  of  the  first  k  parties,  when  the  last  n-k  outcomes  are  un-observed  and  hence  summed  over.  Note,  this  marginal 
probability  distribution  P(al,a2,...,ak  xl,...,xlt )  formally  depends  on  all  n  measurement  settings.  The  last  equality  in 
(5)  imposes  the  NS  constraint  requiring  that  the  marginal  probability  depends  only  upon  the  k  measurement  settings  of 
the  parties  participating  in  the  joint  measurement  (and  not  on  the  remaining  n-k  measurement  setting  of  the  unobserved 
outcomes). 

As  first  pointed  out  by  Pospescu  and  Rohrlich3,  the  NS  constraint  (5)  by  itself  does  not  single  out  classical  and  quantum 
theories,  i.e.  |S|  <  S0.  PR  proposed  the  following  joint  probability  distribution  for  two  parties  (Alice  and  Bob)  with  two 
measurement  settings  (inputs)  x,y  ={0,1},  and  two  measurement  outcomes  (outputs)  a,b  ={0,1}  given  by 


PR  Box:  P(a,b\  x,y)  = 


if  a  ©  b  =  x  ■  y 
otherwise 


(6) 


By  considering  all  possible  inputs  and  outputs,  it  is  straightforward  to  show  that  PR  correlations  of  (6)  satisfy  all  the 
requirements  for  a  NS  theory  as  follows:  normalization  (total  probability) 
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(7) 


X  P(a,b\x,y) 

a,b={  0,1} 

=  m o  I  *, y)  +  P(U  |X, y)  +P{ 0, 1 1  X,y)  +  P(l, 0  |  y, y) , 
a(Bb  =  0  a®b  =  1 

d/2  +  1/2)  ^  +  (1/2  +  1/2)  \x-y  ’ 
=  A,  +  Z>,  , 

0,x-_y  l,X'_r  ’ 

=  1  Vx,y, 

and  the  NS  constraint 


P(a\x,y)=  ^  P{a,b\x,y) 

*={0,1} 

=  P(a,Q\x,y)  +  P(a,l\x,y ) 
a©/>  =  a©0  =  a  a©/>  =  affil  =  a 

=  U2  SatX.y  +  l/2SatX.y,  (8) 

Jl/2  +  0  (if a  =  0& y -  y  =  0),  0  +  l/2(if  a  =  0& Y-y  =  1) 

[0  +  1/2  (if a  =  l&Y-y  =  0),  l/2  +  0(ifa  =  l&Y-y  =  l)  ’ 

=  1/2  Va,x,y, 

=  P(a  |  x)  Va,  y  (  =>  Isotropic,  i.e.  P(a  \  x)  =  1/2  indep  of  a,  y). 

With  the  PR  Box  define  above  in  (6)  we  can  compute  correlations  as 

E(a,b  |  x,y)  =  P( 0, 0  |  x,  y)  +  P(l,  1 1  x,  y)  -P(0,1  \  x, y)  -P(1,0  |  x,  y), 
a®b- 0  a®b= 1 

(1/2  +  1/2)  ^-(1/2  +  1/2 )5u,yt  (9) 

_|+1  if  x-y  =  0,  i.e.(*,y)e{(0,0),  (0,1),  (1,0)}, 

[-1  if  v ■  y  =  1 ,  i.e.(v,y)  =  (1,1), 


where  we  have  used  E(a,b\  x,y)  =  +] ,a\b'  P(a,b  \  x,y),  where  a'  =  l-2a  (//  =  1  —  2Z?)  associates  the 

measurement  values  «'(//)  e  {+1,-1}  with  the  measurement  value  labels  (bits)  a  (h)  e  {0,1},  respectively.  Therefore,  in 
Fig.  1,  assigning  Alice’s  measurement  directions  A,B  =JA  the  bit  labels  xA=l  and  v5=0,  and  Bob’s  measurement 
directions  C,D  =2>  the  bit  labels  yy  =0  and  j’/;=l ,  and  using  (9)  yields  the  algebraic  maximum  5^=4  of  the  CHSH 
inequality,  as  illustrated  in  Fig.  2. 


PR  Box: 

P(a,b\x,y)  = 


if  a  ®  b  =  x  •  y 
otherwise 


with  measurements  x,  y  and  outcomes  a,  b  as  bits , 
i.e  a,b,x,y  e  {0,1}  (Note:  {0,1}  ++{+1,-1}) 


*B-yc=° 

2?  -A 

Oc=0)C  a(BC)  B(.r,=0) 


PR  Box  yield  algebraic  maximum  of  S 

5  =  E(AC)  +  E{BC)+E{BD)  -  E(AD) 
=  1  +  1  +  1  -  (-1) 
=  4 


Fig.2  PR  Box  with  joint  probability  distribution  achieving  the  algebraic  maximum  5^=4  of  the  CHSH  inequality. 
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Since  Sm=4>Sq=2'I2,  no  quantum  (i.e.  physically  realizable)  state  can  reproduce  the  above  PR  probability  (6).  However, 
the  following  ‘‘state”9  O  =  a+  |o+^<t>+ 1  +  a_  |  with  a±  =  (1± \f2)/l  and  Bell  states  1®^  =  (|00)±|ll))/V2 , 

yields  the  PR  probability  (6)  through  the  usual  trace  rule  PPR{a,b  \  x,y)  =  Tr[OMx  ®Mb]  with  {MxaA  ,MX‘}  =  {cr2 ,  cr, }  and 
{M*c  ,  MxbD  }  =  {(cr,  +  cr2  )/V2  ,  (a,  -  cr2  )/-x/2}  ,  where  {cr,  },.E{i  2  3}  are  the  usual  Pauli  matrices.  Note  that  the  form  of  the 
joint  measurement  between  Alice  and  Bob  written  as  a  pure  tensor  product  of  local  observables  Mxa  ®Mb,  ensures  the 
locality  of  the  spacelike  separated  measurements,  which  cannot  increase  entanglement  between  the  parties.  (A 
measurement  involving  the  sum  of  pure  tensor  products,  such  as  Mx  ®  Mb  +  M,x  ®  M'by  which  might  possibly  create 
entanglement,  would  involve  non-local  measurements  between  the  parties,  which  could  only  be  physically  realized  if  the 
parties  were  brought  together).  The  important  point  is  that  O  does  not  represent  a  physical  quantum  state  since  it  is  non¬ 
positive,  i.e.  it  possesses  the  negative  eigenvalue  a_  =  (1- V2)/2  .  Henceforth,  we  shall  refer  to  non-positive,  unit  trace 
Hermitian  operators  O  capable  of  producing  NS  probability  distributions  as  “states,”  and  reserve  the  specific  term 
“quantum  state”  or  “q-state”  for  the  physically  realizable  positive,  unit  trace  Hermitian  operators  denoted  as  p>0,  (i.e. 
density  matrix). 

Following  Acin  et  al.9  we  now  wish  to  investigate  all  sets  of  «-party  spacelike  correlations  in  terms  of  local  quantum 
observables  (measurements)  Mnon _sig  =  M*1  ®---®  Mx"  that  ensure  NS.  These  correlations  can  be  written  in  the  form 


P0=P0(a1,...,a„\xl,...,xn)  =  Tr[OM;i  0---0MJ],  (10) 

Without  loss  of  generality,  we  can  take  the  local  measurement  operators  Mx  =  FT'  =  |a)  (a  |  to  be  the  projection 

operators  onto  “spin-component”  a  in  the  “direction”  x.  Requiring  that  proper  probabilities  be  derived  from  all  local 
quantum  measurements  imposes  the  condition  that  O  be  positive  on  all  product  states.  This  implies  that  0=  W  is  an 
entanglement  witness10  (EW)  with  the  property  0.  Here  some  definition  are  helpful.  A  q-state  is 

separable  (contains  only  classical  correlations)  if  it  is  of  the  form  psep  =  /  ,  p,pf 1  ®pf2  ®--®  pfN  where  each  pf  is  a 
local  density  matrix  and  ^  pi  =  1.  (If  a  q-state  is  not  separable,  it  is  entangled).  Each  local  density  matrix  has  a  (non¬ 
unique)  ensemble  decomposition  pf*  ~  ^  .  pb  |  where  '^Jjpkij  =1-  The  requirement  that  W  is  positive  on  all 

product  states  (a,  |  JR  |  a,  /?,...)>  0  ensures  that  Tr[psepW ]  >  0  from  the  form  of  psep.  A  q-state  p  such  that 

Tr[pW\  <  0  is  then  entangled  (since  it  is  not  separable),  and  JRis  said  to  “witness”  (or  exhibit)  the  entanglement  of  p. 
Note  that  W  is  in  general  a  non-positive  Hermitian  operator.  In  the  context  of  (10),  we  now  consider  0—*W  as  a  state 
(not  necessarily  a  q-state)  from  which  to  derive  NS  correlations  through  the  joint  probability  distributions 

Pw=P(al,...,aN\x1,...,xn)  =  Tr[WM2  ®-®M*:]>  0.  (11) 

The  correlations  (11)  are  termed  Gleason  correlations  by  Acin  et  al.9 

The  subtle  distinction  between  (10)  and  (1 1)  is  that  the  latter  produces  positive  probabilities  for  all  local  NS 
measurements,  while  the  former  may  produce  non-negative  probabilities  on  only  a  subset  of  NS  measurements.  This 
distinction  is  important  since  it  has  been  shown10'1 1  that  for  bipartite  systems  n= 2,  any  Gleason  correlation 
P(a1,a2  \xux2)  =  Tr[W M*1  ®M*2  ]  >  0  can  be  converted  to  a  probability  distribution  derived  from  a  q-state 

=  |  with  modiried  measurements  P(a, ,  a2  \xl,x2)  =  Tr[W  M*1  ®  M*2  ]  =  Tr{p^m)  M*1  ®Mxa;]>  0.  Here 

|  <Dra  'j  is  any  pure  bipartite  state  (not  necessarily  maximally  entangled).  The  proof  relies  on  the  explicit  use  of  the  Choi- 
Jamiolkowski  isomorphism1011’12  (CJI)  which  allows  any  bipartite  ( n=2 )  witness  W  to  be  written  as 
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W("  2)  =  (I®A)(p^^),  where  A  is  a  positive  trace  preserving  map.  In  the  above,  M**  =  A* (Mj) )  where  A*  is  the 
adjoint  of  the  map  A,  i.e.  Tr[A  A(R)]=Tr[A*(H).8].  The  proof  then  follows  directly  as 

Pw  (a,  b  |  x,  y)  =  Tr[W  Ml  0  Myb  ]  =  Tr[(I  0  A)(p  )  Mxa  ®  Ml  ] 

_  (12) 

=  TrWl  ®Mn/®A)(^))]  =  Tr[Ml  ®  A*  (Ml)  p^l  =  Mxa  ®Myb)l 

where  the  second  equality  uses  the  CJI,  the  third  equality  uses  the  cyclic  property  of  the  trace,  the  fourth  inequality 
utilizes  I  ®  A  acting  to  the  left  on  the  tensor  product  measurements  Mxa  ®  Mb  thereby  introducing  the  adjoint  A* 
operation  and  the  modified  local  measurement  operation  Mj)  =  N*  (My )  in  the  last  equality.  Acin  et  al.9  point  out  that 
the  CJI  decomposition  W(N=2}  =  (I  ®  A)(p^^  , )  in  general  fails  for  n> 2  (which  they  demonstrate  by  a  specific  example). 
Thus,  the  Gleason  correlations  (11)  are  strictly  larger  (|.S'|  >Sq)  than  quantum  correlations  for  n>2  (and  equivalent  only 
for  n<  2).  The  state  O  =  a+  |o+^0+ 1  +  a_  |  used  in  the  example  of  PR  correlations  in  the  discussion  after  Fig.  2 

is  not  an  EW  since  it  can  produce  negative  probabilities  for  measurements  other  than  those  considered  (it  would  be  an 
EW  if  it  produced  positive  probabilities  for  all  measurement  choices).  Acin  et  al 9  classify  the  distributions 
P(a],...,an  \xl,...,xn)as  (i)  No-Signaling  if  and  only  iff  can  be  written  in  the  form  of  (10),  (ii)  Quantum  whenever  O  is 
positive  (O>0),  and  (iii)  Local  if  and  only  if  O  corresponds  to  a  separable  quantum  state.  In  the  following,  we  investigate 
the  NS  correlations  of  (10)  and  the  conditions  for  which  they  become  either  Gleason,  or  Quantum  correlations. 

2.3  No  Signaling  (NS)  Correlations:  2-Qubits 

Following  Acin  et  al9  we  define  an  «-partite  probability  distribution  P(al,...,an  \  x1,...,xn)as  being  NS  if  and  only  if 
there  exists  local  quantum  measurements  Mj'  and  a  Hermitian  operator  O  of  unit  trace  such  that  (10)  holds.  It  is 
important  to  note  that  O  need  not  produce  positive  probabilities  for  other  measurements  outside  this  set.  Acin  et  al9  give 
a  prescription  for  the  formal  construction  of  O  given  the  set  of  measurements  Mxa‘ .  In  the  following  we  present  an 

explicit  construction  for  O  for  the  case  of  n= 2  qubits  (r= 2  outputs,  i.e.  a,b  =  {0,1})  and  arbitrary  number  m  of 
measurement  inputs  ( x,y  =  {0,l,...,/«-l }).  Later,  we  extend  this  to  the  case  of  «=3  for  qubits. 

As  stated  in  the  previous  section,  without  loss  of  generality  we  can  take  the  local  Hermitian  measurement  operators  to  be 
the  projection  operators  onto  “spin-component" ^ ’  a  in  the  “direction”  x,  Mxa  =  Tla  =  |  a)  (a | .  For  each  x,  the  completeness 

of  the  measurement  operators  give  y,'  =  Irxr  where  Irxr  =  I  is  the  r  x  r  identity  matrix.  This  allows  us  to  write  the 

a=r-l  measurement  operator  as  =  Irxr  -  ^  ■  One  defines  the  (tilde)  Hermitian  matrices  Mxa  dual  to 

Mxa  through  the  inner  product  Tr\MxaMxa,  ]  =  Sx  x'8a  a,.  For  the  bipartite  case  n=  2,  with  in  general  m  measurement  settings 
with  r  measurement  outcomes,  one  has 


r- 2  m- 1  /■- 2  m- 1  r- 2  m- 1 

0=  X  X  P(a,b\x,y)Ml®Ml+YJYJP(a\x)Ml®I+YJYJP(b\y)I®Ml+I®I,  (13) 

a,b-0  x,y-0  a-0  x=0  b= 0  y 

where  I  is  the  tilde  matrix  dual  to  the  r  x  r  identity  matrix  /,  with  the  additional  orthogonality  conditions  defined  by 
Tr[I I]  =  7>[7]  =  1,  Tr\  M'a  1  \  =  0,  and  Tr\I M\ ]  =  Tr\M\  ]  =  0.  The  conditions  ensure  that  O  is  Hermitian,  Tr[  0\=  1  and 
probabilities  are  given  by  the  trace  formulas  P(a,b\  x,y)  =  Tr[OMxa  0M,f] ,  P(a  \  x)  =  Tr\OMxa  0/]  and  P(b  \  y)  = 

Tr[0 1 0  M l  ] .  This  is  illustrated  in  Fig.  3  where  Alice  and  Bob  share  PR  correlations  by  means  of,  what  are  termed  in 
the  literature,  a  pair  of  PR  boxes  (or  NS  {non-signaling}  boxes). 
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P(a,b  |  x,  v)  =  Tr[OMxa  ®A//]  >  0 


Alice  and  Bob  share 
paired  (generalized) 
PR  Boxes  with 

P{a,b\x,y)  =  \!  2 
itT a®b  =  X' y  mod 2, 
(0  otherwise) 


Fig.3  PR  Box  shared  between  Alice  and  Bob. 


In  the  following  we  specialize  to  the  case  of  qubits  (i~=2,  a,b  =  {0,1})  with  arbitrary  number  m  of  measurement  inputs 
{x,y  =  {0,1,..., 777-1})  .  In  this  case  the  measurement  operators  M*l=0  are  given  as  projection  operators  for  “spin-up”  along 
the  directions  x  —>  mx  on  the  Bloch  sphere.  The  M/=0  are  just  density  matrices  on  the  Bloch  sphere  written  as 


M*=0  =  1 0)v  (0 1  =  1/2  (7  +  mx  ■  a),  mx  =  mx  (sin  Qx  cos  <f>x ,  sin  0X  cos  <j>x ,  cos  0X ) 

\mx  \  <  1,  (density  matrix  on  Bloch  Sphere), 


(14) 


where  a  =  (al,a2,a3)  is  the  vector  of  single  qubit  Pauli  matrices.  Although  not  required  for  the  case  of  qubits,  the 
projection  onto  “spin-down”  along  x  is  given  by  M/=1  =  |l)r  (l|  =  1/2 (7  -mx  -  a)  =  7-M*=0,  with  /  the  2x2  identity 
matrix.  Equation  (13)  now  simplifies  to  the  form 

m-l  m-l  m — 1 

O  =  Y,  p(a  =  0, b  =  0 1  *, y)M 0  ®  Ml  +  Y, P(a  =  0 1  x)Mx0  <S> /  +  ^ P(b  =  0 1  y)I  ®M}0  +I®I.  (15) 

jc,_y=0  x=0  j;=0 

We  simplify  the  notation  by  defining  y,M*=0;x  =  0, -  -  - , m  — 1|  =  { M _j  =  7,{M&0}  =  {A7_, ,  A70,  AT, ,. . .}}  =  {Ma={_l  i>0)}  (a 
set  of  7W+1  linear  independent  matrices)  with  duals  {M/S={_1  y.>0}}  =  |a/_,  =  I,M0,Ml . ..J  satisfying  the  trace  orthogonality 
conditions  Tr[MaM /j]=Sa  p,  and  similarly  for  |/,{Mj'a0}|  — >  {/V>e={_1  We  therefore  write  (15)  as 


m — 1  m — 1  m-\ 

°=  Z  Pij^i  ®^i  +Z^'A4r,  ®I  +  YJP'°I®Nj  +I®I, 

i,j= 0  i=0  7=0 


(16) 


using  the  abbreviations  P°j°  =  P(a  =  0,b  =  0  |  x  -  i,y  =  j) ,  P°’’  =  P(a  =  0 1  x  =  i)  and  P"’°  =  P(b  =  0  |  y  =  y).  For  the 
measurement  matrices  M  x  =  /  =  /2x7,  andM&0  =  \/2(I  +  mi  ■  a),  |i7i(|  <  1,  the  dual  matrices  are  given  explicitly  by 
M  j  =  /  =  1/2(7  >o7M(.  •  a)  =  l/2(7-7w  •  <r),  and  M(.  =  mi  ■  a  where  tt7(.  p  |/«(.  |  >  1,  with  the  orthogonality 

relations  Tr[f]  =  1,  Pr[M .  ]  =  0,  Tr\Mi  7]  =  0,  and  Tr[MM/ 1  =  ®.  Using  the  relationship  Tr\X  ®  Y]  =  7/-[A]  7VfF]  it  is 
straightforward  to  verify  that  Tr[(9]=l  and,  for  example,  P°j°  =Tr[OMi  ®  Nj  ]  which  picks  out  the  term  Mi  ®Nj  in 
(16).  Other  probabilities  are  obtained  for  example  as  P°jl  =  Tr[0  A7(  ®  (7  -  Nj )]  =  Tr[OMi  ®I]-Tr[OMi  ®Nj] 

=  P°’’-P°f  =  Zi,-;nn Pi>f  ~ Pt  f  =  Puj  =  p(a  ~  0, b  =  1 1  x  =  /,  j  =  y)-  Substituting  the  explicit  expressions  for  the  dual 
matrices  into  (16)  yields  the  general  expression  for  O  in  terms  of  products  of  Pauli  matrices 
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(17) 


1  in  - 1 

O  =  -[  Z  (4 C  -  2(7’°’’  +  if0)  + 1)  (m,  •  (?)  0  (h,  •  <?) 

4  f,y=o 

m-1  m-1 

+  ^  (2 if’  - 1)  (m,  •  (?)  0  7  +  J] (2P;°  - 1)  7 0(« .  •(?)  +  /  0 1  ]. 

i=0  7=0 

Specializing  to  the  PR  correlations  in  (6)  given  by  P(a,b  \  x  =  i,y  =  j)  =  1/2  Sa9b  i.j ,mail  =>P°f  =1/2 S()  i  ;mod2  with 
marginals  T}0’*  =  Zf  =1/  2  Vi,  j,  yields  the  expression  for  the  NS  PR  operator 


°pr  =^[  ('?e-C7)0(«e-cr)  +  Owf -O-)0(«o- cr) +  (/«„•  cr)0(«f-cr) -(/«„•  £7)0 («0-<T) +  7  0/  ], 
where  me=  £  m2i,  in0=  £  m2f+1,  we  =  X  &y>  «„  =  Z  tfey+i- 

*=0,1,2,...  i=0,l,2,...  7=0,1, 2,...  y=0, 1,2,. . . 


(18) 


In  (18)  the  subscripts  {e,o}  denote  {even,odd}for  the  summation  over  even  and  odd  dual  measurement  vectors.  Note  that 
in  (18)  the  “single-o”  terms  cr  0/  and  700=  (representing  measurements  by  Alice  or  Bob  alone,  respectively)  have 
dropped  out  since  the  marginal  distributions  P{a\x)=P{b\y)=\/2  are  independent  of  a,b,x,y.  This  leaves  only  the  solely 
two-party  correlation  terms  cr  0cr;  and  the  maximally  mixed  term(7  07)/4.  For  the  bipartite  case  n= 2  often  considered 

in  the  literature  for  two  qubits,  each  with  two  measurement  directions  x  e  for  Alice  and  y  =  {«„,«, }  for  Bob  (i.e. 

a,b,x,y  =  {0,1})  we  obtain  the  simplified  form 

0'PR  =^-[  {m0  -(?)0(«o  •(?)  +  (/«„  •c?)0(«1  •(?)  +  (»;,  -(?)0(«o  •(?)-(/«,  •(?)0(«1  •(?)  +  /  0  7  ].  (19) 


Using  the  procedure  for  calculating  probabilities  discussed  after  equation  (16) ,  the  following  probabilities  can  be 
computed  from  (19) 


«0 

«, 

«0 

«, 

»0 

«, 

»0 

«o  «, 

m0 

1/2 

1/2' 

m0 

'1/2 

1/2' 

m0 

"0 

0  " 

m0 

"0 

0  ' 

m0 

T  l ' 

mx 

.1/2 

0  _ 

mx 

.1/2 

0 

mx 

0 

1/2. 

mx 

0 

1/2. 

mx 

l  -l 

P{a=Q,b=Q\x=mi,y=n ;),  P(a  =1,6  =\\mnnJ),  P(a  =0,b  =l\mj,nJ),  P(a  =  \,b  =  0\mi,nj),  E(m.,nj) 


a®b  =  0, 

■  P  =  l/2for {x, y)  e  {(0,0), (0, 1),(1,0)}, 


a  ®b  =  1, 

■  P  =  l/2for(x,j)e{(l,l)}. 


Here,  the  correlations  in  (20)  are  computed  as  (see  (1)) 

E{m.,n. )(.  .  {0  j,  =  P(a=0,b=0  \  x=m.,y=m. )  +  P{a  =  l,b  =  1 1  iii.Ji.  )  —  P(a  =  0,7>  =  1 1  m.,n )  -  P(a  =  \,b  =  0  j 


(20) 


(21) 


with  corresponding  S  parameter  (see  (2)) 

S  =  E(m0  ,«„)  +  E(m0 , «, )  +  E(mt  ,n0)~  E(fhl  ,n1)  =  4  =  SAM, 


(22) 


achieving  the  algebraic  maximum  value  SiM  =  4. 

For  the  case  of  two  qubits  with  m= 3  measurement  vectors  x  e  {»}„,«;,  m3}  for  Alice  and  j  =  {«„,«,,  n2 }  for  Bob  (i.e.  a,b  = 
{0,1},  with x,y  =  {0,1,2})  we  obtain  from  (18)  the  probabilities  and  correlations 
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(23) 


«o 

>h 

«o 

n2 

no 
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=>  P  =  1/2  for  (x,  y)  e  {(e,  e),  (e,  o),  (o,  e)} ,  ^P  =  1/2  for  (x,  >>)  e  {(o,  o)} . 


In  (23)  e  =  {0,2}  denotes  even  indices  of  the  measurement  directions  while  o  =  {1 }  denotes  odd  indices.  We  achieve  the 
algebraic  maximum  for  the  S  parameter,  generalizing  (22)  defined  as 

S  =  E(me,  ne)  +  E(me  ,na)  +  E(m0 ,«.)-  E(m0  ,n0)  =  4  =  SAM.  (24) 

Note  that  the  dimension  of  the  measurement  vectors  in.  is  set  by  the  dimension  D  =  d2  - 1  of  the  Hilbert  space  of  the 
observer,  which  simply  states  that  any  (D  + 1)  x  (D  + 1)  matrix  can  be  written  in  term  of  the  (D  + 1)  x  (D  + 1)  identity  matrix 
and  the  D  generators  of  su(d).  For  qubits,  D= 3  and  the  three  generators  of  su(2)  are  the  usual  Pauli  matrices  a.  For  a 
given  set  of  m  measurement  3-vectors  \mt\  (vectors  in  the  Bloch  sphere,  |  mi  |<  1 )  one  needs  to  solve  for  the  correspond 
dual  measurement  vectors  {w, }  satisfying  mi  ■mj  =  8i  ..  We  write  these  equations  as  the  matrix  equation 

M,„x3M3xm  =Imxm  where  the  /th  row  (i  =  {0,l,...,w-l})  of  (the  known  coefficient  matrix)  Mmx3  is  mn  and  the  /'th  column 
of  (unknowns)  M3xmis  ilij .  By  linear  algebra,  there  exists  a  right  inverse  of  Mmx3  via  M,is(1Ilj=Mr(MM7’}1  (if 
(MM1’}1  exists)  if  the  columns  of  Mmx3  span  R'n ,  which  can  only  occur  for  m<D= 3.  The  systems  of  equations  is  under¬ 
determined  and  there  exists  at  least  one  solution  (typically  and  infinite  number  due  to  undetermined  free  parameters). 
This  is  the  situation  for  probabilities  and  correlations  shown  in  (20)  and  (23)  for  the  case  m= 2  and  m= 3  measurement 
vectors,  respectively.  For  the  m>D= 3,  there  exists  at  most  one,  unique  solution  (if  any).  This  is  the  least  squares  (LS) 

solution  using  the  pseudo-inverse  Mmx3  given  by  MiS  =  (MrM)_1Mr  (if  (MrM)~‘  exits).  In  general,  the  LS  solution 
has  non-zero  residual  errors  given  by  Err  =  Mmx3  MiS(3xm)  -Imxm,  corresponding  to  joint  probabilities  that  may  be 
negative  for  some  measurements  but  still  satisfy  the  (total  probability)  normalization  condition 
y  t  P(a,  b  |  x,  y)  =  1,  Vx,  y.  Nonetheless,  it  is  instructive  to  perform  numerical  searches  in  the  case  of  m> 3  of  random 

measurement  vectors  to  seek  solutions  which  yield  all  joint  probabilities  in  the  range  Q<P(a,b\x,y)<\,  for  all  pairs  of 
measurement  vectors  mi , n;  for  Alice  and  Bob  that  still  yield  supra-correlations,  i.e.  0  <  S-Sq<  4-2a/2  =  1.172. 

For  the  case  m= 4,  a  particular  solution  is  shown  in  (25)  that  yields  S-Sg=  0.102  (for  brevity,  we  only  show 
P(a=0,b=0  |  x=mi ,  y=n  )  and  the  correlations  E(m.,n.)).  In  general,  the  even/odd  structure  of  the  correlations  E(mj,n.) 
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exhibited  in  the  cases  m< 3  ((20)  and  (23))  is  destroyed,  yet  they  still  produce  supra-correlations  S-So>0.  For  each  value 
of  m  in  Fig.  4  (left)  we  searched  105  random  trials  of  the  measurement  vectors  ,E{0il  and  plot  the  value  of  S- 
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Sq  for  the  first  solution  encountered  in  which  (i)  we  find  proper  joint  probability  distributions  0 <P(a,b  \  x=ni ,  y~nj  )<  1  for 
all  measurement  vectors,  and  (ii)  which  produce  supra-correlations,  S-Sq> 0.  In  Fig.  4  (middle),  we  plot  the  minimum 


Figure  4.  (left)  Numerical  simulations  for  m={2,3,4,...,12}measurement  vectors  yielding  valid  joint  probabilities 
0  <  P{a,b  |  x  =  ifli ,y  =  Hj)<  1  for  all  pairs  of  measurement  vectors  {mj ,  H  - },  /G,0 ,  m  l,  and  supra-correlations 

0.  (middle)  Minimum  eigenvalue  Amln  of  O  in  (18).  The  negative  values  indicates  that  O  is  not  a  proper 
quantum  state  (positive,  Hennitian  operator,  p>0).  (right)  Iteration  number  (out  of  105  random  trials)  at  which  the 
first  set  of  m  measurements  vectors  were  found  for  Alice  and  Bob  producing  supra-correlation  shown  in  the 
leftmost  figure. 

eigenvalue  A.11]in  of  the  matrix  O  in  (18).  The  negative  value  of  /.rnm  indicates  that  O  is  not  realized  by  a  proper 
quantum  state  (i.e.  a  positive,  Hermitian  operator,  p>0).  The  rightmost  plot  in  Fig.  4  is  the  iteration  number  at 
which  the  first  set  of  measurement  vectors  was  found  which  produced  supra-correlations.  For  the  values  of  13  < 
m  <  20  numerically  explored,  no  supra-correlations  solutions  were  found  within  10s  trials  (the  plot  indicates  that 
it  becomes  exponentially  hard  to  find  such  a  solution). 

2.4  No  Signaling  (NS)  Correlations:  3-Qubits 

The  bipartite  results  of  the  previous  section  for  «=2-qubits  are  straightforwardly  extended  to  the  tripartite  case  of  n= 3- 
qubits  with  similar  implications.  Flere  the  generalization  of  the  bipartite  CFISFl  nonlocality  parameter  S  is  given  by  the 
SvetlichnyL’  inequality  (SI)  relating  correlations  E(a,b,c\x,y,c)  between  three  spacelike  separated  parties  A,  B,  C 

S  =E(a,b,c  |  0,0,0 )  +  E(a,b,c  |  0,1,0 )  +  E(a,b,c  \  1,0,0) -E(a,b,c  |  1,1,0) 

+  E(a,b,c  |  0,0,1) -E(a,b,c  |  0,1,1) -E{a,b,c  \  1,0,1)- E(a,b,c  \  1,1,1).  ^ 

The  SI  has  the  bounds  (i)  |S|<SC  =  4  for  classical  correlations,  (ii)  |i>|  <SQ=  4x[2  for  quantum  correlations,  with  (iii) 

the  algebraic  upper  bound  given  by  |S|<^m=8,  achieved  when  the  correlations  in  (26)  take  the  values  E=  1  if  they  are 

preceded  by  a  plus  sign,  and  E=- 1  if  they  are  preceded  by  a  minus  sign.  The  generalization  of  the  PR  correlations  of  (6) 
is  given  by14 


TPR  Box:  P(a,b,  c\x,y, z) 


Jl/4  if a®b®c  =  x-y®y-z®x-z 
[  0  otherwise 


(27) 


often  referred  to  as  a  tripartite  PR  (TPR)  box.  The  marginal  distributions  of  (27)  are  again  isotropic  and  satisfy  the  NS 
constraint,  i.e.  P(a,b\x,y,z)  =  P(a,b\x,y)=H4  for  all  a,b,x,y,z  and  P(a,  \x,y)=P(a\x)=H2  for  all  a,x,y,  and  similarly  for  all 
other  marginal  probability  distributions. 

For  the  case  of  n= 3  qubits  (r=  2  output  measurement  values)  a,b,c  ={0,1),  with  m  possible  measurement  vectors  for  each 
observer,  x,y,z  ={0,l,...,m-l)  we  again  find  that  only  the  highest  (three  party)  correlations  term  and  the  maximally 
mixed  term  are  non-zero  in  the  expression  for  Ojpr 


Proc.  of  SPIE  Vol.  8400  84000Y-10 
10 

Downloaded  from  SPIE  Digital  Library  on  15  May  2012  to  132.3.29.68.  Terms  of  Use:  http://spiedl.org/terms 


Otpr  =  U{{me-d)®{ne-d)  +  {me-d)®{n0-d)  +  (m0  ■  a)  ®  («,  •  (?)  -  (m„  ■  (?)  0  (n0  •  cr)}  0  (re  •  a) 

O 

-  {(m„  •  (?)  ®  (n0  •  (?)  +  (m0  •  (?)  0  («(,  •  (?)  +  (me  •  (?)  0  («0  •  (?)  -  (m(,  •  (?)  0  (ne  •  (?)}  0  (r0  ■  (?)  +  /  0 1 0  /],  (28) 

where  qe  =  £  q2i,  q0=  Y  4+ i>  |  =  {w,  n,  r}. 

~  *=0,1,2,...  ~  ~  i =0,1,2,...  ~ 

The  regular,  even/odd  (mod  2)  structure  of  0TPR  in  (28)  reflects  the  non-zero  structure  of  the  TRP  probabilities  in  (27), 
and  can  be  seen  as  an  additional  single  qubit  generalization  of  0PR  in  (18).  That  is,  the  2-qubit  term  in  the  first  curly 
brackets  in  (28)  {(ine  -d)  0  (ne  •(?)  +  ( me  -d)  0  (na  ■<?)  +  ( mo  -d)®{ne  ■d)-(mo  -d)®{n0  •(?)}  tensor-producted  with  the 
remaining  “even”  qubit  term  (re  •  (?) ,  is  precisely  two-party  correlation  term  that  appears  in  0PR  in  (18).  Similarly,  the 
term  in  the  second  curly  bracket  in  (28)  {(mo  ■  d)  0  (na  ■  d)  +  (mo  ■  d)  0  (ne  ■  (?)  +  (ine  ■  (?)  0  (no  ■  d)  -  (me  ■  (?)  0  (ne  ■  d)} 
tensor-producted  with  the  remaining  “odd”  qubit  term  (7o  ■ d )  (with  the  accompanying  minus  sign)  is  just  the  bit  flip 

(e<->o)  of  the  previous  two-party  correlation  term.  Again,  we  can  achieve  the  algebraic  maximum  ,S'  )  W=8  when  each  party 
has  (for  the  case  of  qubits)  at  most  m= 3  measurement  vectors  (for  exactly  the  same  linear  algebraic  reason  for  the  n= 2 
bipartite  case).  Further,  as  in  the  bipartite  case,  we  can  find  particular  NS  supra-correlation  solutions  0<5-5o<4-2a/2  for 
m> 3,  but  which  become  increasingly  hard  to  find  the  larger  the  value  of  m. 

3.  POSITIVITY  AND  PURITY  CONSTRAINTS  ON  QUANTUM  CORRELATIONS 

in  general,  any  «-quibt  2"  x  2”  matrix  can  be  decomposed  into  the  sum  of  tensor  products  of  at  most  n  Pauli  matrices  with 
the  2x2  identity  matrix  I.  For  the  bipartite  case  composed  of  n= 2  qubits  we  can  write  the  matrix  Ou,‘2)  as15,16 


0(2)  = 


4 


I  ®  I  +  (sm 


•  (?)  0  /  +  /  0  (?(2)  •  (?)  +  Y  Cap  <y a  ®  c p 

a,p= 1 


(29) 


Fleres01  and  s(2)  are  the  vectors  on  the  Bloch  sphere  representing  Alice’s  and  Bob’s  single-party  properties  respectively, 
and  the  3  x  3  matrix  Ca/j  =  Tr[0 oa  0  a p  ]  encode  the  2-party,  joint  correlations  between  Alice  and  Bob.  Flere  we  have 
used 


=SMvT  +  is^(J^  and  TrUl 2]  =  1.  (30) 

From  (17)  and  (18)  s0>  =J(2)=0  with  the  PR  correlation  matrix  given  by  =  ^'”.'n(4P°:°  -  X)(mi)a  ®(nj)p,  where 
P°j°  =  1/2  §0  ,  and  hence  (477°:°  —  1)  =  1  if  /  -  j  mod  2  =  0,  and  - 1  otherwise.  The  essential  feature  of  PR  correlation 

matrix  is  reflected  in  the  coefficient  (4 P°j°  -1)  involving  the  PR  probabilities,  which  reproduces  the  rightmost 
3x3  matrix  E(m. ,  n. )  in  (23)  for  the  case  m= 3  (which  is  the  well-posed  case  with  3  equations  for  3  unknowns  given  the 
mj ,  for  the  dual  measurement  vectors  in ,  satisfying  mt  ■  ini  =  A  . ,  possessing  a  unique  solution).  Thus,  it  is  reasonable  to 
drop  the  requirement  that  the  set  of  fixed  measurement  vectors  {w/ }  and  {«,  }  are  predetermined,  and  agreed  upon  by  Alice 
and  Bob,  and  consider  instead  the  solely  two-party  correlated  state 


i ®i+  Y  cpJva®°, 
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fxPR  _ 
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-l  l 

,  =^>  E(m,n)  =  ^  maC^  ®  np  =  m  -CPR  -n. 

(31) 

a,fi= 1 

l 

l 

l 

a, 

Equation  (31)  yields  supra-correlations  that  achieve  the  algebraic  maximum  SAM= 4  with  the  choice  e.g.  {mi={0 , 2, }  = 
{x,y,zj  =  {«  ={01 2}}  (amongst  many  possibilities).  Flowever,  for  arbitrary  measurement  vectors  {»;,.}  and  \nt }  the 
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correlations  are  bounded  by  |£’(w,«)|  <  2V2  (in  numerical  simulations,  by  2.551  in  105  random  trials),  which  can  lead  to 
negative  probabilities  since  e.g.  P(a  =  0,7  =  0  |  x=m,y=n)  -  1/4(1  +  E{m,n)).  If  one  scales  C™  -»  C™  /  (2^2)  so  that 
\E(m,  n)\  <  1,  the  probabilities  are  now  bounded  between  0  and  Vi,  and  0PR  becomes  an  entanglement  witness  (a  non¬ 
positive  Hermitian  operator  that  is  positive  on  all  product  states),  but  an  uninteresting  one,  since  the  value  of  the 
nonlocality  parameter  is  reduced  from  Sam=4  to  SAM  12^2  =^2<SCl=2. 


As  relevant  comparison  to  0PR  ,  the  density  matrix  for  the  quantum  singlet  Bell  state 
by 


3 

i 

1 

O 

O 

7®7+  Z  C^aa®cjp 

singlet  _  c 

9  ap  ~  ~~  °ap  ~ 

0-10 

Ct,p=\ 

0  0-1 

^singtet^  _  (|oi)-|lo))/V2  is  given 
=>  Esmgla  (in,  n)  =  -in  ■  n,  (32) 


(note,  the  other  three  maximally  entangled  Bell  states  also  have  diagonal  2-party  correlation  matrices  with  two  +1,  and 
one  -1  matrix  elements).  yosmglet  is  a  valid  quantum  (pure)  state  with  eigenvalues  {1,0, 0,0},  while  0PR  is  not,  since  it 

possess  a  negative  eigenvalue  (1  -  ~J\1  )/4  .  The  question  at  hand  is  what  governs  the  structure  of  the  2-party  correlation 
matrix  Cu„  ?  In  the  following,  we  consider  the  constraints  of  both  positivity  (i.e.  non-negative  eigenvalues)  and  purity 
(i.e.  is  quantum  state  pure)  on  0PR. 


For  completeness,  we  note  that  the  tripartite  PR  box  Otpr  in  (28)  is  of  the  solely  three-party  correlation  form 

l[ 


0TpR 


7®/®/+  X  CZoa®a,®o 

a,p,y=\ 

which  is  a  special  case  of  the  n= 3  qubit  state  0<,!=3) 


1 


(33) 


=  —  [  7  ®  7  ®  7  +  ( J  •  <j)  ®  /  ®  /  +  /  ®  .  <j)  ®  7  +  7  ®  7  ®  (s''  ’  •  (?) 

8 

+  Z  +  Z  BaP°a®I®°p  +  Z  CafiI®CTa®°e+  Z  C afiy  01 G P  0  1' ■ 

a,p-\  «,/?=!  cc,p- 1  a,p,y-\ 


(34) 


Again,  the  pertinent  question  is  what  governs  the  distribution  of  the  correlations  amongst  single  parties  (  s(,) ),  two 
parties  Aafj ,  Bafj ,  Cap ,  and  three-parties  CaPjt . 


3.1  Positivity  constraints 

In  an  extensive  review  of  2-qubit  states  of  the  form  (29),  Englert  and  Metwally 16  derive  constraints  amongst  the  single 
and  two-party  correlations  imposed  by  the  (positivity)  condition  that  0(2)  >  0  .  Here  we  generalize  their  arguments  to  an 
arbitrary  n-qubit  state  O1"1  defined  by 


1  3 

=_  [/®/®...®/  +  (j11>.(t)®/®...®/  +  ...  +  /®.--/®(J(")-(t)+  Z  cra®ap®I®---®I  + 

N  a,,a2=  1 

-+  z  7®-®/®^,®^+...+  z 

an-Xa„  =1  axa2,..;<xn=  1  ^ 


(35) 


where  IN  =  7  ®  7  ®  •••  ®  7  is  the  Ax N  identity  matrix  with  A  =  2",  and  K  is  the  traceless  part  of Oln>  (noting  Tr\o ,]  =  0) 
given  by 
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K=I-NO(n)  (N  =  2n),  0(n>  >0=>K<1. 


(36) 


K  contains  all  the  terms  in  O1"’ except  the  maximally  mixed  term  I N  /  N .  If  0M  >  0  (all  eigenvalues  are  non-negative) 
then  the  eigenvalues  of  K  are  all  less  than  unity  by  (36).  Now,  the  eigenvalues  {Ai=i  w}  of  any  NxN  complex  matrix  A 

are  determined  by  its  characteristic  equation  Za(A  =  ]~[ ^  (A  —  At )  =  0.  By  the  Cayley-Hamilton  theorem17,  A  satisfies 
its  own  (matrix)  characteristic  equation,  i.e.  Za(^)  =  0-  This  can  be  used  to  develop  a  set  of  recursion  relations11'  relating 
the  coefficients  of  the  characteristic  (Mh  degree)  polynomial  /4(/l)to  traces  of  the  powers  of  A  denoted  by  Tm  =  Tr[A"‘] 


ZAW-r+^r '+...+ dn_xa  +dn 


inDm  +  DmlTt  +  Dm 


-2T2  +  ■ 


.  +  D,Tm_i+Tm  =0,  m  =  \,2,  —  ,N.  (37) 


The  last  term  in  (37)  allows  one  to  recursively  solve  for  Dm  as  rnDm  =  TtDm  k  (defining  D0  =  0).  Specializing  to 
the  traceless  case  A=K  of  (36),  we  have  7j  =  Tr\K\  =  0.  For  the  «=2 -qubit  case  (1V=4)  we  have 

n  =  2,  N  =  2"  =4:  D2=-\/2T2,  D3  =  -1/377.,  D4  =  -l/4[T4  -1/277,2].  (38) 


Note  that  expressions  Dm  for  an  arbitrary  number  of  qubits  n  do  not  require  a  recalculation  of  previously  computed 
expression  of  the  Dm  for  fewer  number  ri  <  n  of  qubits. 


We  now  note  the  following  useful  observation:  since  zA  ( K )  =  0,  and  by  (36)  each  eigenvalue  is  less  than  unity 
{Ai=i  N  <  1},  it  follows  that  ZAW  =  [  [ .  (A  — A(.)  >  0  and  is  monotonically  increasing  for  A  >  l.This  further  implies  the 
same  is  true  for  all  derivatives  of  the  characteristic  polynomial  i.e.  (d  ZaWI dAk)m  >  0,  for  k  =  0,1,. ..,1V.  Thus,  taking 
derivatives  of  the  characteristic  polynomial  written  in  the  form  ZAW  =  ZT-o  DN_kAk  and  evaluating  at  A  -  1  allows  us  to 
define  N  non-trivial  positivity  conditions  (for  derivatives  k  =  2,...,N)  in  the  compact  form 


z 


"N-A 

yN-kj 


(~D)< 


N\ 

k\(N-k)\’ 


(39) 


to  determine  when Oin)  >  0  for  the  case  of  «-qubits.  Using  (39)  with  the  relations  in  (38)  for  £=2,3,4  yields 


T2< 


N(N  - 1);  ( N  -2 )(T2 /2)  +  (T3 /3)  <  T  j ;  (N  ^  3)  (T2 /2)  +  (N- -  3) (T3 /3)  + 1/4 (T4  - 1/2 T2 )  *  f  J • 


(40) 


The  D2  constraint  on  7)  =  7>[  K 1  \  (first  inequality  in  (40)),  can  be  further  simplified.  For  the  general  «-qubit  state  O1"'  the 
non-zero  contributions  to  T2  come  from  the  square  of  each  individual  term,  yielding  a  common  factor  of  N  =  7>[7V  ] 


T2  =  Tr[K2]  =  N 


■  n  =  2:  s2W+s2(2> 


S  (cs?)  +...+  I 

j=1  i< j  ax,a2=\  ax  a2,...,an  = 1 

+  Z  C2<(A7-1W=3, 


<N(N- 1) 


(41) 


a,  ,a-f  =1 


3 


3 


■n  =  3:  J2,1)+?2(2)+J2(3)  +  Y  A2  +  Y  B2  +  V  C2  +  V  C2  <(N-l)  ,=7, 

Au  axa2  A-^  a\ai  A-^  a\a2  A—i  «i«2a3  ^  ' N=2 3  ’ 

a1,a2=  1  ax,a2  =  1  ax,a2  =  1  ax,a2,a3=l 
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where  we  have  illustrated  the  specific  cases  of  n= 2  from  (29),  and  n= 3  from  (34).  Equation  (4 1 )  can  be  intuitively 
understood  as  follows.  The  first  line  of  (4 1 )  indicates  that  there  are  N  — 1  =  2”  - 1  ways  to  distribute  correlations  in  the 


state  O' 


i («) . 


=  n  single-party  correlations  s  (1),s  t2)  ...,s 


=  n(n  — 1)/2  two-party  correlations  for  !</</<« 


etc...,  up  to 


=  1  /7-body  correlations.  The  sum  of  all  these  correlation  repositories  is  y//, 


r  n\ 


ykj 


=  2“  -\  =  N 


We  can  now  make  some  quantitative  observations.  From  (31),  the  n= 2  qubits  PR  state  0{2)  has  s (1)  =  s  (2)  =  0,  with  all 
the  correlations  in  the  two-party  correlation  term  C™  containing  9  non-zero  entries.  Thus,  (c™ )  =  9  violates  the 

n= 2  T2  positivity  constraint  y3  (c„„)  <  3  in  (41).  The  supra-correlation  PR  joint  probability  distribution  in  (6)  gives 
rise  to  more  correlations  in  C™  that  are  allowed  physically  by  nature.  On  the  other  hand,  note  that  the  quantum  singlet 
state  in  (32)  with  C/^glet  =  -8afj  not  only  satisfies,  by  also  saturates  the  T2  positivity  constraint  (as  do  the  other  three 
maximally  entangled  Bell  states). 

As  a  counter  example,  for  n= 3  the  state  0TPR.e  from  (33),  defined  by  setting  =  £aPr  (the  Levi-Civita  anti-symmetric 

symbol  with  value  +1  for  even  permutations  of  (1,2,3),  and  -1  for  odd  permutations)  has  6  non-zero  entries,  and 
therefore  does  satisfies  the  «=3  73  constraint  y  3  ,(C^;e )  =6<7  in  (41).  However.  0TPRe  is  a  non-positive 

operator,  and  thus  not  a  quantum  state.  To  discriminate  the  non-positivity  of  this  state,  we  must  turn  to  the  Z)3  positivity 
constraint  (the  second  inequality  in  (40);  k=  3  in  (39)).  In  computing  7)  =  Tr[K2\  for  Otpr  one  must  form  the  product  of 
three  Pauli  matrices  for  each  qubit  o^cr^oy  =  isaPrI  +  8a/ja„  -  8ayap  +  SpytJu  with  Tr[aa<jpcrr ]  =  2 isaPy, ,  which  is  purely 
imaginary.  For  the  odd  number  n= 3  qubits  in  OIVRl: .  the  product  of  three  such  traces  in  the  calculation  of  T,  is 

proportional  to  i 3  =  —i .  Since  K  is  Hermitian  with  real  eigenvalues,  the  same  is  true  for  all  positive  powers  K"' ,  so  that 
Tm  in  general  must  be  real.  Thus,  symmetry  relations  in  the  triple  summations  involved  in  computing  7)  conspire  to  yield 
T3  =  0.  The  same  is  true  for  the  trace  of  any  odd  power  of  A  in  0TPR.e,  i.e.  T2m+1  =  7r[7f2"'+1]  =  0,  since  (i)  for  each  qubit, 
traces  of  odd  number  of  Pauli  o-matrices  are  proportional  to  i  and  (ii)  for  an  odd  number  2n’  + 1  of  qubits,  the  product  of 
an  odd  number  of  such  individual  qubit  traces  is  proportional  to  i2n  +1  which  is  again  imaginary.  In  fact,  T2m+l  =  0  is  true 
for  any  odd  number  qubit  state  O12" +l1  that  involves  solely  (2/7'  + 1)  -party  correlations  (i.e.  only  the  first  and  last  terms  in 
(35)).  Substituting  T}  =  0  into  (40)  leads  to  the  more  restrictive  constraint  T2  <  N(N-l)/3  which  translates  into 

Si  p  \{Capr)  -  7/3,  which  is  violated  by  ^  ,(C™;s)  =  6  for  0TPRe  .  Since  the  tripartite  PR  state  is  of  the  form 
0TPR  in  (33)  one  might  try  to  achieve  the  algebraic  maximum  of  the  Svetlichny  inequality  by  requiring  all  8  terms  in  (26) 
to  be  of  equal  magnitude  c.  The  above  constraint  would  require  (Caj8j,)  =  8c2  <  7/3  or  c  =  1/2  -^7/6  if  the  bound 

is  saturated.  But  then  the  Svetlichny  inequality  would  be  bounded  by  Sa  <\S  \=  8c  =  4^7/6  <  S0.  Thus  the  correlations 
could  be  quantum,  i.e.  stronger  than  classical,  but  not  supra-correlations  (stronger  than  quantum). 

3.2  Purity  constraints 

The  77-qubit  state  O1”'  (35)  is  called  pure  if  it  is  a  rank  one  projector,  i.e.  0R'Jre  =  |'PA,//vPiV  |  for  normalized  A-vector 
|'PJV)  ,  or  equivalently,  the  condition  02Ju"l  =  O/"/  holds,  i.e.  ~IN)  =  0.  Upon  writing  OtH)  =  1/N[In - K ]  the 

purity  condition  becomes 
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Pure  State:  0(£(0£-IN)=  0  =>  K2  =  (N-1)In  -{N-2)K. 


(42) 


T  aking  the  trace  in  (42 )  yields  T2  =  N(N  -  1)  which  saturates  the  upper  bound  of  the  first  (D2)  positivity  inequality  in 
(40).  Multiplying  (42)  by  K"'  and  taking  the  trace  yields  the  recursion  relation  Tm+2  =(N-  Y)Tm  -(N-  2)TmAl  which  can 
be  used  to  express  Tm  as  functions  of N\  T}  =  -N(N -T)(N -2),  TA  =  N(N-Y)2  +  N(N -V)(N -2)2 ,  etc...  These  in  turn 
can  be  used  to  explicitly  show  that  the  pure  states  saturate  the  upper  bound  of  the  positivity  constraints  (40),  and 
consequently  all  the  positivity  constraints  (39).  As  an  immediate  consequence,  states  such  as  Otpr  in  (33)  and  their 
generalization,  involving  an  odd  number  n  =  2 «'  + 1  of  qubits  and  solely  «-party  correlations,  cannot  be  pure  due  to  the 
fact  that  all  odd  traces  T2n,+l  =  0  from  the  discussion  in  the  previous  section.  This  is  not  necessarily  the  case  for  states 
involving  solely  n  -party  correlations  when  n  is  even.  Important  particular  instances  are  the  «=2-qubit  maximally 
entangled  Bell  states  involving  solely  2-party  correlations,  of  the  which  the  singlet  state  in  (32)  is  a  representative 
example. 

The  purity  constraint  (42)  forces  relationships  amongst  the  correlations,  determining  how  they  are  distributed  amongst 
the  parties.  An  explicit  construction  of  (42)  for  the  general  2-qubit  state  Oa>  (29)  yields  the  relations  s ^  Caps(p , 

42>  =Zl=i5«)C^’and  C«P  =s«>s(/)- Cafib) where  CT  s  (the  matrix  of  cofactors  of 

Cap  )•  For  the  Bell  states  s®  =  s'p  =  0,  so  that  the  purity  condition  requires  the  2-party  correlations  to  obey 
Gap  =  -Cpb) .  Englert  and  Metwally16  have  shown  that  the  generic  n=2 -qubit  quantum  state  0(2)  =  pa)  >  0  (29)  is  of 
the  form 

n  =  2  qubit  pure  quantum  state:  ppre  =  \/4[I®I  +  p<Jl®I  +  pI®ai-al®al-qa2®a2-ai®a3]>0,  (43) 

where  is  the  common  length  of  the  individual  Block  vectors  s  and  s  <2\  and  q  =  p\  -  p2  >  0  is  the  concurrence19 

(a  complete  measure  of  two  qubit  entanglement  for  both  pure  and  mixed  quantum  states).  For  the  case  of  pure  quantum 
states  (43)  the  concurrence  is  given  by  q2  =  Tr[p(p'ep(2Jre],  where  p<2>'(,  is  the  complex  conjugate  of  p{pre. 

For  n=3  qubit  pure  states  0{ 3 ’  with  general  from  (34),  the  distribution  of  the  correlations  amongst  the  various  parties  is 
much  more  complex.  From  the  discussion  above  we  know  the  correlations  in  (34)  cannot  be  solely  3-party  correlations, 
implying  that  some  of  the  correlations  must  also  be  distributed  amongst  1 -party  and/or  2-party  combinations.  Recently, 
much  research  has  concentrated  on  the  study  of  the  relation  between  nonlocality,  embodied  in  the  Svetlichny  parameter 
S  in  (26)),  and  pure  state  entanglement  as  measured  by  Wootters  tangle20  r.  The  tangle  between  qubits  1,  2  and  3  is 
defined  by 


0  <  r  =  02m)  -  £22  -  &22  <  I .  (44) 

In  (44)  the  concurrence19  quantifies  the  bipartite  entanglement  between  qubit  1  and  2  after  qubit  3  has  been  traced 
out,  and  similarly  for  ^l  v  The  concurrence  Sm2)  measures  the  bipartite  entanglement  between  the  two  subsystems 
composed  of  (i)  qubit  1  and  (ii)  qubits  2  and  3  taken  together.  The  reduced  density  matrix  p2i  =  Tr\pm]  can  be  shown  to 
have  two  non-zero  eigenvalues  and  hence  acts  as  an  effective  single  qubit  whose  bipartite  entanglement  with  qubit  1  is 
measured  by  Sll23)  =  2N/7>[  /?,  |  (for  the  case  of  a  pure  3-qubit  quantum  state).  The  concurrence  £  between  two  qubits 
described  by  the  quantum  state  p  is  given  by19  £  =  max{0,/t,  -A2-A3-A4}  where  X  are  the  square  roots  of  the 
eigenvalues,  in  decreasing  order,  of  the  matrix  pp  where  p  -  (cr2  ®  o2)p\o2  ®a2). 

Ghose  et  al.21  have  investigated  the  relationship  between  the  maximum  value  Smax  of  the  Svetlichny  parameter  for  a 
given  value  of  the  tangle  x  for  a  class  of  three  parameter  3-qubit  states  defined  by  \<f)  =  cos  9X 1 000)  +  sin  9 , 1 )  with 
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|$)  =  |l),  \<j>2)  =  cos  ft  |0/  +  sin#2  |l/,  and  |^3)  =  cos#,  |0/  +  sin#3  |l)  .  The  state  |^)  has  tangle  r  =  (sin  2#,  sin  ft  sin#,)2. 
Some  important  special  cases  of  these  three  parameter  states  are  (i)  the  maximally  sliced  (MS)  states,  of  which  there  are 
two  types,  (MS',)  and  (MS',)  and  (ii)  the  generalized  GHZ  state  |  GGHZ /  given  by 

|MS2)  =  -j=(|000)  +  cos#2 1 1 0 1)  +  sin #2 1 1 1 1)) ,  \MS3)  =  -j=(  |000)  +  cos#3  |l  10}  +  sin#,  |ll l)); 

| GGHZ)  =  cos#,  1 000^  +  sin  #,  1 1 1 1),  with  t2ms  =  sin2  #2,  t2ms^  -  sin2  #,,  t2gghz  -  sin2  2#,. 


Emary  and  Beenaker22  have  established  bounds  on  Smax  for  a  given  three  tangle  as  |  —  5,2ax  - 1 1<  t  <  for  the  MS 

16  32 

states.  Ghose  et  al.21  have  established  the  bound  as  =  4v/l  +  r(MS) .  For  the  GGHZ  state  Ghose  et  al.2'  have 
established  the  bounds  S£™Z)  =  4v/l-r<G®z,forr  <  1/3  and  S(™m  =  4v,2r(GG/iZ)forr  >  1/3.  All  three  states  in  (45) 
achieve  the  maximum  tangle  x=l  on  the  GHZ  state  |  GHZ /  =  l/  V2  (|  000^  + 1 1 1 l/)  (the  generalization  of  the  symmetric 
Bell  state  j®  )  =  l/V2  (|00)  +  |ll)) )  which  also  achieves  the  maximum  nonlocality  S=Sq=4^2  for  a  quantum  state. 


To  examine  how  the  correlations  are  distributed  amongst  the  parties  for  the  states  in  (45)  let  us  denote  sf  =  s2'1'  for  the 
square  of  the  1-party  correlations  in  (34),  Cj2  =  A =  Tr[ArA\,  C,23  =  Tr[BTB] ,  C;3  =  Tr[CT C]  as  the  sum  of  the 

squares  of  the  2-party  correlations,  and  C123  =  C2^  as  the  sum  of  the  squares  of  the  3-party  correlations.  The 

state  | MS',  /  has  correlations  distributed  such  that  (i)  sf  =  sf  =  0,  sf  =  1  -  r,  (ii)  Cj2  =3-2 r,  Cf3  =  C]3  =  r,  and  (iii) 

C23  =  3  +  r,  with  t  -  sin2  #,.  Maximum  tangle  t  =  1  on  the  GHZ  state  arranges  the  correlations  so  that  (i)  all  the  1 -party 
correlations  are  zero,  (ii)  C12  =  C23  -  C22,  =  1  and  (iii)  C,2,,  =  4  is  maximized.  This  is  the  same  behavior  as  the  state 
'Zsvetiicnv }  =  l/ 2(—  1 000^  + 1 0 1 1/  + 1 1 0 1/  + 1 1 10))  that  Svetlichny13  put  forth  as  an  example  of  a  state  that  achieves  maximum 
nonlocality  S=S0 .  This  is  in  contrast  to  another  important  class  of  3-qubits  states,  the  W  state23,  with  r  =  0,  defined  by 
W  =  l/V3(|l00/  +  |010^  +  |001^)  in  which  (i)  the  only  non-zero  1-party  correlations  are  s31)  =  s,32)  =  s\3)  =  1/3,  (ii)  all  2- 
party  correlation  matrices  are  identical  C12  =  C„  =  C23  =  diagonal[2/3,2/3,-l/3]  and  (iii)  the  non-zero  3-party 
correlations  are  given  by  (C123)  =  2/3  for  afiy  =  { 113,131,223,232,312,321}  and  {Cm  ^  = -1.  For  the  generalized 

GHZ  state  \GGHZ}(i)  s,31)  =  s32)  =  s’,3'  =  Vl - r  with  t  =  sin2  2#,  so  that  s,2  +  s;  +^32  =  3 -3r,  (ii)  C12  =  C13  =  C23  = 
diagonal[0,0,l],  with  C,2  +  C23  +  C23  =  3,  and  (iii)  C123  =  1  +  3r.  We  note  that  the  cusp  in  S{°°HZ)  that  occurs  at  r  -  1/3 
discussed  above  occurs  at  the  condition  sf  +  sf  +  sf  =  C123 .  For  the  state  | MS3 )  this  latter  condition  would  yield  the 

unphysical  result  r  =  -1 .  The  lesson  to  be  learned  from  this  non-exhaustive  examination  of  3-qubit  pure  states  is  that 
besides  classifying  tripartite  entanglement  solely  in  terms  of  the  tangle,  there  is  discriminating  information  to  be  gleaned 
from  the  distribution  of  the  correlations  amongst  1-,  2-,  and  3 -parties. 

4.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  have  examined  the  structure  of  supra-correlations  that  are  stronger  than  quantum  and  hence  not 
realizable  by  a  physical  (positive)  quantum  state  p>0.  The  supra-correlations  are  intriguing  because  they  arise  from  valid 
probability  distributions,  first  put  forth  by  Popescu  and  Rohrlich  (PR),  that  satisfy  the  no-signaling  principle  of  special 
relativity  as  well  as  all  the  usual  normalization  condition  on  the  joint  and  marginal  distributions.  Thus,  the  fact  that 
nature  is  not  able  to  realize  these  supra-correlations  points  to  hidden  structure  underlying  how  quantum  correlations  can 
be  distributed  amongst  spacelike  separated  parties.  This  paper  has  examined  the  structure  and  distribution  of  PR 
correlations  in  2-  and  3 -qubit  systems  by  explicitly  constructing  “states”  (not  necessarily  positive  quantum  states)  that 
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exhibit  supra-correlations  for  a  fixed,  but  arbitrary  number,  of  measurements  available  to  each  party.  We  have  shown 
that  the  PR  correlations  involve  only  solely  n -party  correlations  amongst  the  n  observers.  We  have  extended  this  study  to 
include  n-party  correlations  that  capture  the  essential  features  of  the  PR  correlations  and  do  not  rely  on  predetermined 
measurements  between  the  n  participants.  By  constructing  constraints  based  on  the  positivity  and  purity  of  an  arbitrary 
n-qubit  state  we  have  shown  the  “unreasonableness”  of  the  PR  correlations  in  that  they  encode  more  correlations  than 
are  physically  allowed  by  nature.  In  future  work  we  will  couple  this  approach  of  studying  how  correlations  are 
distributed  amongst  the  n  parties  to  the  study  of  quantum  entanglement.  The  study  of  entanglement24  is  an  important,  but 
difficult  field,  only  well  understood  for  the  case  of  two  qubits  (both  pure  and  mixed),  and  to  a  lesser  degree,  for  pure  3- 
qubit  systems.  A  fruitful  area  to  investigate  next  are  pure  3 -qubit  systems,  where  a  generalized  (though  non-unique) 
Schmidt  decomposition  holds25.  We  purport  that  an  examination  of  the  distribution  of  correlations,  bounded  by 
physically  imposed  constraints  on  e.g.  positivity  and  purity,  coupled  with  the  description  of  entanglement  in  terms  of  the 
tangle,  as  initiated  in  this  work,  can  shed  further  light  on  the  classification  of  pure  tripartite  systems. 
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